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1 Basic Statistical Mechanics 

This introduces the basic ideas using a system of distinguishable particles, 
e.g. magnetic ions in a paramagnetic salt. The initial goal is to find the heat 
capacity C. We can get this if we know the total energy U . 

Each particle has energy levels e±, e 2 , ... with ni, n 2 , ... particles respec- 
tively. 

In real systems, there may be constraints like total particle number N 
and energy U must be fixed. 

For given constraints, we assume that the particles are equally likely to 
have any arrangement (or microstate) among the energy levels. 

We want to find 

U = ni£i + n 2 e 2 + 

so we need to know the most likely energy distribution (rii, n 2 , ...) (or macrostate). 
"Most likely" means that number of arrangement 



ni\n 2 \... 

is largest. We need to find the unknown {rii,n 2 , ...) when Q is maximum. 

This is tricky when we have constraints on N and U. We need to use the 
Lagrange multiplier method. Instead of Q, we find the maximum of 

L = lnfi + AiiV + \ 2 U 

where Ai, A2 are additional unknowns. 



1 



Using Stirling's approximation, 



In n! ~ n In n — n 



we then differentiate L with respect to rii and solve for rij. For just this 
differentiation, we assume that terms N and U in L are not fixed and are 
functions of rii (but iV in the Q function is still fixed). 
Together with the Boltzmann postulate 



An important implication is that at the microscopic level, temperature can 
be determined by particle populations (ni, n 2 , ...) - e.g. negative temperature 
(e.g. laser), magnetic cooling. 

To get U, we also need E.g. in a spin-1/2 salt, these are ~HbB and 
+HbB. Substituting these into Boltzmann distribution and solving for A 
gives 



In an ideal gas, we assume that the particles do not interact with each other. 
Because of the large number of energy levels £j, it may be convenient if the 
sum for U can be approximated to an integral. So we need a way to find the 
number of states in de. 

We solve Schrodinger's equation for a particle in a cubic box of side a. 
The solution can have the form 



S = k B InQ 



and the entropy relation 




the final solution is Boltzmann distribution 





2 Ideal Gas 



ip = A sin k x x sin k y y sin k z z 
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We need only positive k x ,k y , k z because the negative values only changes the 
sign of the wavefunction and gives the same state. The other parts of the 
solution are 

h 2 k 2 

e = — — and k 2 = kl + k 2 + k 2 

2m x y z 

Assuming that potential is infinity at the walls, wavefunction goes to zero 
there. This results in the quantisation conditions 

i, _ Ux7T i. _ n y n Ah- nz71 
a a a 

The possible points in k space give a 3D grid of points it /a apart. If we divide 
the space into identical cubes centered at every point, we can associate the 
volume (tt/o) 3 with each point. 

By considering the volume of a spherical shell of radius k and thickness 
dk, we can count the number of states in the shell: 

g(k)dk = 47rk 2 dk -T- (vr/a) 3 + 8 

This defines g(k), the density of states, k is related to e. In the corresponding 
de, there are the same number of states: 

g(e)de = g(k)dk 

This gives the density of states in energy variable 

. AirmV f- 

g{e) = h3 V2me 

If we know the number of particles f(e) in a state, the total energy is 

roc 

U= / ef(e)g(e)de 



f(e) corresponds to n*. We can find it using the Lagrange multiplier method 
again. This time, the particles are indistinguishable. 

We consider arrangements of particles and states in de, obtain an expres- 
sion for all energy and maximise this with constraints on iV and U. For a 
typical gas at room temperature, the number of particles in any de is much 
smaller than the number of states. This assumption leads to the Boltzmann 
distribution again: 

/(E )=Aexp(-^ 
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The integral above can now be solved to give 



U = ^N A k B T 



for 1 mol of the gas and heat capacity 

C = ^N A k B . 

3 Electrons in Metals 

A simple model of electrons in metals assumes that the conduction electrons 
form a gas. It neglects the rest of the atoms. This model looks like an ideal 
gas, but the measured electronic heat capacity of a metal is much smaller. 

The difference is because of Pauli's exclusion principle. At room temper- 
ature, electrons in metals stack up in energy levels. There is not free states 
above for them to get excited to - unless they are very near the top - the 
Fermi level. 

To find C, we can follow similar steps to ideal gas. The density of states 
is doubled because there are 2 spin states for each energy state: 

. , n AirmV r- — 
g(e) = 2 x — — — V2me 
h 6 

We find the number of arrangements of particles among states in each interval 
de following similar steps as ideal gas. But we do not assume that the number 
of particles in any de is much smaller than the number of states. We must 
also insist that no 2 particles can be in the same state. 

Then applying the Lagrange multiplier method with the constraints on 
N and U gives the Fermi-Dirac distribution: 

m 



exp[(e - fi)/k B T] + 1 

where \x is an unknown constant that depends on N and T. At T = 0, f(e) 
is 1 when e is below fi and 0 otherwise. Then number of states below /i must 
equal number of particles. The relation can be written as 



! g(e)de = N 
Jo 
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Then we can solve for \i. Its value at 0 K is called Fermi energy: 



f37T 2 N\ 

F 2m 



To find the heat capacity, we estimate the number n of electrons excited 
above E F . Obtain the approximate distribution above E F 

/(£) = e ^BT e -e/k B T 

The width of k B T means that 

n w g(E F )k B T 

Since this /(e) has the same form as Boltzmann distribution, the ideal gas 
result applies. So total energy of excited electrons is 

C/i = n x 

Combining these results gives the molar heat capacity as 

dU _ $k B T 
dT 2 Ep 



4 Phonons and Photons 



4.1 Phonons 

Heat capacity of a solid also comes from vibration of atoms. This can be 
modelled as sinusoidal waves of displacements of atoms from rest positions. 

If we assume that these displacements are zero at the sides of a solid 
cube, we get the same quantisation conditions as ideal gas. This mean we 
can use the same density of states. But we must multiply that by 3 because 
for each wavevector (k x , k y , k z ), there can be 2 transverse and 1 longitudinal 
polarisations. 

Another difference is that angular frequency uj is often used instead of e. 
The relation with wavevector is 



co = v k 
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where v is speed of sound in the solid. Transform density of states formula 
using 

g(oj)duj = g(k)dk 

to get 

Vcu 2 

^ )=3x w 

Waves with frequency uj are quantised into packets of energy 

e = hw 

Think of packets as particles - called phonons. Obtain their energy dis- 
tribution f(e) using similar steps as for electron gas but without exclusion 
principle. Also no constraint on particle number because phonons are virtual 
particles. Lagrange multiplier method then gives 

^ = exp{hu/k B T) - 1 

No n because no constraint on N. Integrate huf(u)g(u) to find U . But 
upper limit is not infinity - because cannot have wavelength smaller than 
spacing between atoms. Expression for total energy is 

run 



U — I fkjf(u)g(uj)du! 
Jo 



for some maximum ojd called Debye frequency. To obtain this, use maths 
result that number of possible modes for iV atoms is 3iV. (A mode is a wave 
with a particular polarisation and wavevector.) Then 



3N = I g{u)duj 



Solving gives 



67VttV\ 1/3 



V 



Combining the above gives 

JJ 3Vh f- CD *:*(!*: 



2 7r 2 t ,3 j o exp(hufkBT) — 1 

At the high and low temperature limits, this can be solved. The resulting 
heat capacities are respectively 

C = 3Nk B and C = AT 3 

where A is an expression containing the constant factors in the above integral. 
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4.2 Photons 



A hole on a hollow box acts like black body in absorbing all radiation falling 
on it. The energy interacts with the inner walls, achieve thermal equilibrium 
and gets remitted through the hole. This is blackbody radiation. Instead of 
heat capacity, we are interested here in its frequency spectrum. 

The steps are very similar to phonons. Instead of atomic vibration, we 
have electromagnetic waves here. Assuming a perfectly conducting box, elec- 
tric field parallel to inner walls must be zero. Assuming sinusoidal waves, 
this boundary condition gives the same quantisation conditions as ideal gas. 

This should mean the same density of states (DOS). But because the wave 
of each wavevector can have 2 transverse polarisation, the ideal gas DOS must 
be multiplied by 2. Like phonons, we work in variable u. Transforming to 
this gives 



where v is now speed of light c. Wave packet energy and f(oo) are the 
same as phonons. But instead of finding the total U, we want the frequency 
distribution of energy u{uj). This the energy per unit frequency: 



Liquid helium is used for cooling by evaporating it. Lowest temperature is 
limited to about 1 K - its vapour pressure falls exponentially and it cannot 
evaporate much below 1 K. 

A mixture of helium-3/helium-4 can cool down to mK. Below about 1 K, 
the mixture separates into 2 layers. The phase more concentrated in helium-3 
about the phase more dilute. Cooling happens when helium-3 diffuses from 
upper to lower layer. 

This downward movement corresponds to the upward movement of vapour 
when liquid helium evaporates. The concentration of helium-3 in lower layer 
corresponds to vapour pressure of liquid helium. Below a few mK, the sat- 
uration concentration of lower layer remains fixed at 6.6%. So helium-3 can 




u(u) = huf(u)g(uj) 



Vhoo 3 1 
7r 2 c 3 exp(hufkBT) — 1 



which is the blackbody radiation spectrum. 



5 Dilution Cooling 
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continue to diffused downwards as long as the lower layer is not saturated. 
This is the advantage over evaporating liquid helium whose vapour pressure 
quickly goes to zero below 1 K. 

Two ways to explain why cooling takes place: 

1. Increase in entropy in lower layer. 

First, obtain the heat capacity for each layer. Helium-3 atoms are fermions 
like electrons. Use the formula from electron gas for molar heat capacity: 

7T 2 T 

C=——R 
2 T F 

where Fermi temperature T F is Ep/ks- Compute heat capacities for each 
layer using correct effective mass ratios and molar volumes.. 

Integrate dS = dQ/T to find entropy. Use 3rd law of thermodynamics to 
set the lower limit, and use the C formula to integrate. Result is 

S = C 

This shows that molar entropy Sd of lower layer is higher than that of upper 
layer Sc- So when helium-3 flows to lower layer, increase in entropy requires 
heat absorbed from surrounding - which then cools down. 

Combining the above, heat change when 1 mol of helium-3 flows to lower 
layer is 

AH = T(S D - S c ) = 84T 2 J/K 

2. Increase in excited particles. 

Since helium-3 are fermions, we can use the expression for number of 
excited electrons 

3iV 

n « g(E F )k B T = ——k B T 

Lower layer has lower Ep because of lower concentration N/V. So it has 
more excited helium-3. 

Average energy of excited helium-3 is the same 3/csT/2 in both layers. So 
total energy of excited helium-3 is higher in lower layer. In moving to lower 
layer, helium-3 must gain energy. It has to absorb this from the surrounding 
- which then cools down. 
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6 Magnetic Cooling 



This method can cool down to microKelvin. Consider a paramagnetic salt 
first. The possible steps are: I. Start at liquid helium temperature; 2. 
Magnetise isothermally; 3. Demagnetise adiabatically. 

Step 2: Consider a spin- 1/2 salt. Increasing applied field increases spac- 
ing between magnetic energy levels. As higher level moves up, particles fall 
to lower level to obey Boltzmann distribution. So heat is given out. Mag- 
netising it slowly gives the heat time to dissipate so temperature stays the 
same. 

Step 3: Then thermally isolate the salt. When field decreases, energy 
spacing decreases. As higher level moves down, particles should move back 
up from lower level. But they cannot - because they cannot get energy to do 
this since the salt is thermally isolated. According to 



if ?7-2 is fixed, then B/T is fixed. So as field decreases, temperature must fall. 

If B falls to zero, T cannot go to zero because of field from neighbouring 
magnetic ions. This field determines the Curie temperature T c below which 
the ions align themselves and become ferromagnetic. Then the salt stops 
working because it cannot respond to the external field. So T c gives a lower 
limit for the cooling. 

For paramagnetic salt, T c is in the milliKelvins. So it is not possible to 
reach microKelvin using this salt. Another material whose ions give much 
smaller field is needed. The nucleus would do. The nuclear magneton is 
about 2000 times smaller than the Bohr magneton for electron. So the field 
would be 2000 times smaller also. Then T c would be well below microKelvin. 
Typically, a block of copper is used in nuclear cooling. 

In practice, it takes a long time - days, weeks - to reach the the final 
temperature. To speed things up: 

1. Nuclear cooling starts at milliKelvin, e.g. using a dilution refrigeration. 

2. It uses much more power magnet. The small \x n also mean small energy 
spacing. Stronger field to increase this spacing would reduce the particles at 
the higher level. When demagnetised, fewer particles at higher level means 
lower temperature (for the same final field). 
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7 Superfluid 

8 Superconductor 
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